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Abstract 

We apply Tsallis's g-indexed nonextensive entropy to formulate a random matrix 
theory (RMT), which may be suitable for systems with mixed regular-chaotic dynam- 
ics. We consider the super-extensive regime oi q < 1. We obtain analytical expressions 
for the level-spacing distributions, which are strictly valid for 2x2 random-matrix en- 
sembles, as usually done in the standard RMT. We compare the results with spacing 
distributions, numerically calculated for random matrix ensembles describing a har- 
monic oscillator perturbed by Gaussian orthogonal and unitary ensembles. 



^email : alaajibdelhady@yahoo.com 
■^email : ayamagd@hotmail.com 



1 INTRODUCTION 



The past decade has witnessed a growing interest in Tsahis' non-extensive generahzation 
of statistical mechanics [1, 2]. The formahsm has been successfully applied to a wide class 
of phenomena; for a review see, e.g. [2, 3]. The standard statistical mechanics is based 
on the Shannon entropy measure S = — J dxf{x)ln f{x) (wc use Boltzmann's constant 
ks = 1), where /(x) is a probability density function. This entropy is extensive. For a 
composite system A + B, constituted of two independent subsystems A and B such that the 
probability p(A + 5) =p{A)p{B), the entropy of the total system S{A + B) = S{A) + S{B). 
Tsallis proposed a non-extensive generalization: Sg = (1 — / dx[f{x)]'^)/ {q — 1). The 
standard statistical mechanics is recovered when the entropic index q = 1- The value of q 
characterizes the degree of extensivity of the system. The entropy of the composite system 
A + B, the Tsallis' measure verifies 

SgiA + B) = Sg{A) + Sg{B) + (1 - q)Sg{A)SgiB), (1) 

from which the denunciation non-cxtensivc comes. Therefore, Sq{A + B) < Sq{A) + Sq{B) 
if q > 1. This case is called sub-extensive. U q < 1, the system is in the super-extensive 
regime. The relation between the parameter q and the underlying microscopic dynamics is 
not fully understood yet. 

Non-extensive generalization to the random matrix theory (RMT) [4] has recently re- 
ceived considerable attention [5, 6]. The conventional RMT is a statistical theory of random 
matrices H whose entries fluctuate as independent Gaussian random variables. The matrix- 
element distribution is obtained by extremizing Shannon's entropy subject to the constraint 
of normalization and existence of the expectation value of Tr(H"'"H) [7]. Non-extensive gen- 
eralizations of RMT, on the other hand, extremize Tsallis' non-extensive entropy, rather 
than Shannon's. The first attempt in this direction is probably due to Evans and Michael 
[8]. Toscano et al. [9] constructed non-Gaussian ensemble by minimizing Tsallis' entropy 
and obtained expressions for the level densities and spacing distributions. A slightly differ- 
ent application of non-extensive statistical mechanics to RMT is due to Nobre and Souza 
[10]. Recently, Bertuola et al. [11] generated a new family of ensembles that unifies some 
important applications of RMT. In the range -oo < q < 1, this family was found to be a 
restricted trace ensemble that interpolates between the bounded trace ensemble (Chapter 
19 of [4]) at q = — oo and the Wigner Gaussian ensemble at q = 1. The non-extensive 
formalism was applied in Ref. [12] to ensembles of 2 x 2 matrices. Analytical expressions 
for the level-spacing distributions of mixed systems belonging to the orthogonal-symmetry 
universality class are obtained in [12]. The calculation of the spacing distribution showed 
different behavior depending on whether q is above or below 1. It is found that the sub- 
extensive regime of g > 1 [13] correspond to the evolution of a mixed system from chaos 
modelled by the standard RMT to order described by the Poisson statistics. On the other 
hand, the spectrum behaves in a different way in the super-extensive regime, where q < 1. 
Starting with a value of g = 1, where the nearest-neighbor spacing (NNS) distribution 
is well approximated by the Wigner surmise (see below), and decreasing q the distribu- 
tion becomes narrower, and more sharply peaked at spacing s = 1. It develops towards the 
picked-fence type, such as the one obtained by Berry and Tabor [14] for the two-dimensional 
harmonic oscillator with non-commensurate frequencies. 

The present paper considers the departure from chaos in the super-extensive regime. 
Section 2 reviews the non-extensive formulation of RMT. It contains the derivation of the 



level spacing distributions for super-extensive systems with and without time reversible sym- 
metry. Section 3 is devoted to a numerical experiment in which a two-dimensional harmonic 
oscillator is perturbed by a random-matrix ensemble. It is shown that the super-extensive 
RMT provides a reasonable description of the final stage of the stochastic transition. The 
conclusion of this work is given in Section 4. 

2 Non-extensive RMT 

RMT replaces the Hamiltonian of the system by an ensemble of Hamiltonians whose matrix 
elements are independent random variables. Dyson [15] showed that there are three generic 
ensembles of random matrices, defined in terms of the symmetry properties of the Hamil- 
tonian. Time-reversal-invariant quantum system are represented by a Gaussian orthogonal 
ensemble (GOE) of random matrices when the system has rotational symmetry and by 
a Gaussian symplectic ensemble (GSE) otherwise. Chaotic systems without time reversal 
invariance are represented by the Gaussian unitary ensemble (GUE). The dimension /3 of 
the underlying parameter space is used to label these three ensembles: for GOE, GUE and 
GSE, j5 takes the values 1, 2 and 4, respectively. Balian [7] derived the weight functions 
Pa(H) for the three Gaussian ensembles from the maximum entropy principle postulating 
the existence of a second moment of the Hamiltonian. He applied the conventional Shannon 
definition for the entropy to ensembles of random matrices as 5 = — /dHP^(H) InPg(H) 
and maximized it under the constraints of normalization and fixed mean value of ^^(H+H) 
and obtained 

P^ir], H) = Z^^ exp [-tjTV (H+H)] , (2) 

which is a Gaussian distribution with inverse variance I/277. Here Z^^ is a normalization 
constant. 

The non-extensive RMT applies the maximum entropy principle, with Tsallis' entropy, 
to obtain matrix-element distributions that are no more independent. The Tsallis entropy 
is defined for the joint matrix-element probability density P^(g, H) by 

S,[Pp{q,m= (1-/ dH[P^(g,H)]^y(g-l). (3) 

We shall refer to the corresponding ensembles as the Tsallis orthogonal ensemble (TsOE), 
the Tsallis Unitary ensemble (TsUE), and the Tsallis symplectic ensemble (TsSE). For 
q 1, Sq tends to Shannon's entropy, which yields the canonical Gaussian orthogonal, 
unitary or symplectic ensembles (GOE, GUE, GSE) [4, 7]. 

There are more than one formulation of non-extensive statistics which mainly differ in 
the definition of the averaging. Some of them are discussed in [16]. We apply the most 
recent formulation [17]. The probability distribution P^(g, H) is obtained by maximizing 
the entropy under two conditions, 



dHP^(g,H) 

JdH [P^(g,HrTr(H+H) 
/dH[P^(g,Hr 



= 



(4) 
(5) 



where the constant is related to the variances of the matrix elements. The optimization 
of Sq with these constraints yields a power-law type for Pp{q,H), which may be expressed 
as [12] 



P^iq, H) = [l - (1 - q)rjg {tV (H+H) - a}} 



1/(1-9) 



(6) 



where r]q = r]/ / dH H)]"^, rj is the lagrange multiplier associated with the constraint 



in Eq. (5), and Zg = / dH [l - (1 - q)?ig |Tr (H+H) - aj^ 
stant. We use in Eq. (6) the notation [u]-^- = max{0, u}. 



1/(1-9) . 



is a normalization con- 



+ 



We now calculate the joint probability density for the eigenvalues of the Hamiltonian 
H. With H = U~''-XU, where U is a unitary matrix. For this purpose, we introduce the 
elements of the diagonal matrix of eigenvalues X = diag(xi, • • • , xn) of the eigenvalues and 
the independent elements of U as new variables. Then the volume element (3) has the form 



dH = \ An (X)|^dXd/x(U), 



(7) 



where A at (X) = Un>mix n Xfm^ is the Vandermonde determinant and (i/i(U) the invariant 
Haar measure of the unitary group [4]. In terms of the new variables, Tr(H"'"H) = Y^f:=i so 
that the right-hand side of Eq. (7) is independent of the angular variables in U. Integrating 
Eq. (6) over /x(U) yields the joint probability density of eigenvalues in the form 



P/3{q,xi,-- ■,xn) = Cp{q) 



W {Xn Xjji) 



n>m 



N 



(8) 



here C/j is a normalization constant. 

Analytical expressions for the nearest-neighbor spacing (NNS) distribution can be ob- 
tained for the 2x2 matrix ensembles. For the Gaussian ensembles, this approach leads to 
the well-known Wigner surmises 



P(s{s) = a/3/exp (-6/35^) , 



(9) 



where {a^,bp) = (f,f), (f, ^) and {^,^) for GOE, GUE, and GSE, respectively. 
They present accurate approximation to the exact results for the case of iV ^ 00. 

We consider the non-extensive generalization of the Wigner surmises, hoping that the 
results present a reasonable approximation to the physically interesting cases of large N, as 
they do well in the case of Gaussian ensembles. We rewrite Eq. (8) for the case of TV = 2, 
and introduce the new variable s = \xi — X2\ and X = (xi + X2)/2 to obtain 



Pp{q,s,X)=Ci:s{q)s^ 



1 - (1 - (l)Vg { + 2^' 



1 

1-q 



(10) 



Now, we consider the case oi q < 1. For this case, the distribution in Eq. (9) has to be 
complemented by the auxiliary condition that the quantity inside the square bracket has to 

be positive. Thus, we integrate over X from —z to z, where z = y^l — ^(1 — q)rjqs'^. We 
obtain the following expression for NNS distribution in the super-extensive regime 



Pf3{q,s) = a(s{q)s^ l-hp{q)s 



1-q 



(11) 



where 



and 



Ml) 



(/3+l)/2 



r(2 + f + ^ 



r(^)r(i + ,^) 



r 



/3+2 



r(2 + f + ^ 



_r(m)r(| + f + ^) 

Thus, in the case of conserved time reversal symmetry when /3 = 1, 



(12) 



(13) 



h (g) 



+ 3Ui(g), 



1-q 

r ("5 + ^ 



In the absence of time reversal symmetry, P = 2, and 

4r (3 + 



(14) 
(15) 



02(9) 

^^2 (g) 



1-9 



[^'2(g)] 



3/2 



4 

TT 



(16) 
(17) 



3 Numerical experiment 



It has been suggested in [12] that the super-extensive regime corresponds to the evolution 
of the chaotic system towards a regular regime like that of a two-dimensional irrational 
oscillator described by Berry and Tabor [14]. The energy eigenvalues for such a system is 
given by 

Ei = hut (n + am) , (18) 



where uj is the larger frequency and a is the frequency ratio while i stands for the quantum- 
number pair m and n. Berry and Tabor calculated 10 000 eigenvalues and constructed the 
histograms of P{s) for oscillators with a = l/\/2, 1/^5, and e~^. These histograms are 
compared in Fig. 1 with the spacing distribution Pi(0,s), calculated using Eq. (11). 
The figure shows that the distribution Pi(0, s) does not exactly fit the numerical histograms. 
This situation here is similar to the case of sub-extensive regime considered in [12, 13]. It 
was found in these papers that the non-extensive RMT cannot reproduce the Poisson 
distribution that described generic regular system for any value of g > 1. It can only 
describe the evolution of the shape spacing distribution until its peak reaches roughly half 
the distance from the peak position of the Wigner distribution to that of the Poisson. 

Wc shall demonstrate that the super-extensive distribution describes the final stage of a 
transition from a regular harmonic dynamics to chaos in a similar way as the sub-extensive 
theory describes the final stage of the stochastic transition of generic regular systems [12, 13]. 

To show this, we shall now consider the following matrix Hamiltonian 



H(5) = (1 - 5)Hho + 5HGOE, 



(19) 




Figure 1: Comparison between NNS distributions for two-dimensional harmonic oscilla- 
tors with different frequency ratios a obtained by Berry and Tabor [14], and the NNS 
distribution of the superstatistical ensemble with g = 0. 

where Hho = diag(£'j) is a diagonal matrix whose diagonal elements are given by Eq. (18) 
while Hgoe is a random matrix of equal dimension with entries drawn from a GOE. By 
varying the parameter g from to 1, the Hamiltonian H((7) describes a transition from a 
regular harmonic-oscillator type regime to chaotic dynamics. 

We have considered ensembles of 20 matrices of dimension 200 x 200 representing the 
Hamiltonian in Eq. (19). We numerically diagonalized these matrices and calculated the 
NNS distributions for different values of the parameter g. We compared the resulting 
spacing distributions to Eq. (11) for /3 = 1 and evaluated the values of the entropic index 
parameter q that best fit these distributions. The results are shown in Fig. 2. We see 
that for g ~ 0.02 we approach the GOE distribution. Fig. 3 shows the variation of the 
values of the entropic index q that best fit these distributions with the strength of the GOE 
perturbation g stating from 5 = where g = and the spacing distributions take the form 
of picked-fence type up to values of g where q approaches 1 and the spacing distributions 
become that of GOE. 

We have considered perturbing the two-dimensional oscillator with a GUE, where the 
Hamiltonian takes the form 

H(5) = (1 - 5)Hho + sHgue, (20) 

Fig. 4 shows the NNS distributions of the two-dimensional oscillator perturbed by a GUE 
of strength g. We compare the resulting spacing distributions to Eq. (11) for (3 = 2 and 
evaluated the values of the entropic index q that best fit these distributions. Fig. 5 shows 
the variation of the entropic index q with the GUE perturbation strength g. 

From Fig. 2 and Fig. 4, we see that Eq. (11) does not exactly fit the initial stage of 
the transition from the regular two-dimensional oscillator dynamics but fits the final stage 
of the transition into both GOE and GUE. From Fig. 3 and Fig. 5, we see that the GUE 
limit is approached faster than the GOE limit. 



4 Conclusion 



In summary, the NNS distribution obtained from applying Tsallis statistics to RMT maps 
two routes for the transition from chaos to order. One leads towards the integrability 
described by the Poisson statistics by increasing g > 1, but ends at the edge of chaos. The 
other is followed by decreasing q form 1 to 0. It leads to a picket-fence type spectrum, such 
as the one obtained by Berry and Tabor [14] for the two-dimensional harmonic oscillator 
with non-commensurate frequencies. 
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Figure 2: NNS distributions for the Hamiltonian matrix representing a two dimensional 
harmonic oscillator perturbed by a GOE Eq. (19) for different values of the parameter g. 
Dashed curves represent fitting the histograms to Eq. (11) for /3 = 1. The values of the 
entropic index q that best fit these distributions are given. The solid curves are pure GOE. 
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Figure 3: Variation the entropic index q of Eq.(ll) for /3 = 1 that best fits the NNS 
distributions of the Hamiltonian in Eq. (19) with the parameter g. 
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